Ab initio calculations of structural and electronic properties of CdTe clusters 
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We present results of a study of small stoichiometric CdnTen (l<w<6) clusters and few medium 
sized non-stoichiometric CdmTen [{m,n = 13, 16, 19); (m^n)] clusters using the Density Func- 
tional formalism and projector augmented wave method within the generalized gradient approxi- 
mation. Structural properties viz. geometry, bond length, symmetry and electronic properties like 
HOMO-LUMO gap, binding energy, ionization potential and nature of bonding etc. have been 
analyzed. Medium sized non-stoichiometric clusters were considered as fragments of the bulk with 
Tfj symmetry. It was observed that upon relaxation, the symmetry changes for the Cd-rich clusters 
whereas the Te-rich clusters retain their symmetry. The Cd-rich clusters develop a HOMO-LUMO 
gap due to relaxation whereas there is no change in the HOMO-LUMO gap of the Te-rich clusters. 
Thus, the symmetry of a cluster seems to be an important factor in determining the HOMO-LUMO 
gap. 

PACS numbers: 61.46.-w, 61.46.Df, 71.15.Mb, 73.22.-f 



I. INTRODUCTION 

Semiconductor nanoparticles or Quantum Dots (QDs), 
in particular of II- VI materials, have received tremendous 
attention during the last two decades owing to their un- 
usual physical properties and wide range of applications-'^. 
A systematic study of QDs is useful to understand the 
evolution of their physical and chemical properties with 
size. Moving from a molecule to bulk, one can observe 
a range of variation of fundamental properties. This size 
dependent effect, however, is most significant for small 
sized clusters and is attributed to the different geometry 
of clusters in comparision to the bulk. Quantum confine- 
ment of electronic states and large surface to volume ratio 
modify the properties of these clusters. In bare clusters, 
unsaturated bonds lead to peculiar structural and elec- 
tronic properties like surface reconstruction, formation of 
new cleavage planes etc. 

The differences in the structural properties also affect 
the electronic properties of the QDs. In bulk systems 
one observes quasi-continuous electronic levels forming 
bands whereas in clusters the electronic levels are dis- 
crete. Apart from the structural and electronic proper- 
ties, other properties like the thermodynamic and optical 
properties too differ significantly for clusters as compared 
to bulkS*^. 

Among the II - VI semiconductors, CdX (X = S, Se, 
Te) have potential application in optical and optoelec- 
tronic devices. CdTe nanoclusters, in particular, are used 
in optical sensors. Multilayered CdTe with poly (dial- 
lyl dimethyl ammonium chloride) are electro and photo- 
active and hence are used in neuroprosthetic devices^. 
Self assembled CdTe nanodots have high efhciency in 
terms of photoluminescence and hence are very good for 
the purpose of making LEDs. CdTe and Au nanohy- 
brid materials have photoluminescence variation depend- 
ing upon environmental condition like temperature, ionic 
strength, solvents and hence are used as biosensors^. 

We present Density Functional Theory (DFT) based 



calculations for some small, Cd„Te„ (l<n<6), and few 
medium sized CdmTen {m,n = 13,16,19) (m^n) clus- 
ters. Both Cd and Te in bulk phase exist in hexag- 
onal structure with c/a ratio 1.89 and 1.33 respec- 
tively^ whereas CdTe in bulk phase has lattice consatnt 
a = 6.48A in zincblende structure and a = 4.57A and 
c — 7.47A in wurtzite structure. CdTe has a calculated 
direct band gap of 1.59eV at the F pointX. 

Study of small clusters is challenging due to lack of ex- 
perimental structural information and increased degree 
of freedom. Very small clusters (containing up to 8-10 
atoms) exhibit symmetry and have only a few possible 
geometries. On the other hand, medium-sized clusters 
can be considered as fragments of bulk possessing the 
same structural symmetry as that of the bulk^. To the 
best of our knowledge, so far no ab-initio calculations 
have been reported on CdTe QDs. We compare our re- 
sults for small CdTe QDs with those of CdS and CdSe 
wherever available. 

Structural, electronic and optical properties of II- VI 
semiconductor nanoclusters have been studied using a 
variety of methods based on ab-initio DFT, quantum 
chemistry methods like Hartree- Fock, Gaussian etc. and 
approximate methods like tight-binding metho d^ ^^ . 



II. METHODOLOGY AND COMPUTATIONAL 
DETAILS 

The first step for calculating the cluster properties is 
the determination of the lowest energy configuration. Ex- 
istence of several configurations with varying coordina- 
tion and multiple minima in the potential energy surface 
enormously conrplicate the problem making it almost im- 
possible to decide by direct calculation of the lowest en- 
ergy configuration for clusters. 

In this work we have used Born-Oppenhcimer molecu- 
lar dynamics within the Kohn-Sham density functional 
framework-''^. The electronic structure is calculated 
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self-consistently using projector augmented wave (PAW) 
methodic as implemented in VASP packageii^ within 
the framework of GGA instead oiVanderhilt ultra-soft 
pseudo-potential(US-PP)i& as it provides better first 
principle simulations as compared to US-PP^^. In pseudo 
potential methods, the effect of core electrons and nuclei 
is replaced by an effective ionic potential and only the 
valence electrons which are directly involved in chemical 
bonding are considered. 

Various parameterizations have been used for different 
exchange-correlation potentials'^ in the DFT formalism 
in the literature. It has been reported that generalized 
gradient approximation (GGA) gives much better results 
for binding energy and band gap as compared to local 
density approximation (LDA) and therefore we have used 
GGA, as developed by Pcrdew-Berkc-Ernzerhof (PBE)i& 
for our calculation. PBE is known to improve total en- 
ergy and molecule atomization energy and therefore the 
electronic properties, as compared to PW9l'^ and the 
results agree with that of the experiments very well. 

The geometry of small clusters varies significantly from 
that of the bulk. We have chosen several possible ini- 
tial geometries by exploiting the symmetry considera- 
tions and by exchanging the positions of anionic and 
cationic atoms. The number of initial structures is small 
for very small clusters. But for clusters with 10-12 atoms, 
lot many configurations have been used to determine the 
lowest energy structure. 

For structure optimization we have used a super-cell 
of size a = 25A and the plane wave cut-off was set to be 
274.3 eV. These values are chosen after performing cal- 
culations for different cell sizes and cutoff energies. We 
have employed conjugate gradient (CG) technique im- 
plemented via Kosugi algorithm-*^ (special Davison block 
iteration scheme) for optimizing geometry. 

The valence electron configurations used for Cd and 
Te are 5s^4d'° and bs^'ap^ respectively. The Ad levels in 
Te atom are well separated from the 5 s level and hence 
can be included in the core. 

CdTe in bulk has zincblende (ZB) structure with 
T;; symmetry. Therefore for medium sized non- 
stoichiometric clusters Cd„iTe„ {m=/=n), the initial geom- 
etry was considered as a fragment of the ZB structure. 
The size of the super-cell was chosen to be a = 30 A which 
is sufficiently large to minimize inter-cluster interactions 
even for the largest cluster size we consider in this work. 
In all the structure minimizations, the force and energy 
convergence obtained was ~ 10~^ eV/A and ^ 10^'' eV 
respectively. 



III. RESULTS AND DISCUSSIONS 

A. Structural properties of Cd„Te,i 

The geometries for lowest energy (LE) of Cd„Te„ ob- 
tained by CG are shown in Fig. ^^-nd the first local min- 
ima (FLM) in Fig. [3 It is interesting to note that for all 
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FIG. 1: Lowest energy configuration of CdnTe„ for (l<n<6) 
The grey atoms are Cd and the black atoms are Te. 



the calculated lowest energy structures, the chalcogenide 
atoms are at the peripheral positions. This is a direct 
consequence of Coulomb repulsion between the lone-pairs 
situated on the chalcogenide atoms which is minimized 
by separating the chalcogenide atoms towards the periph- 
eral position. As expected, the ground state structures 
of these CdTe clusters do not resemble the bulk phase. 

CdTe dimer is a linear molecule with bond length of 
2.57A and belongs to Coov point group. For Cd2Te2, 
CdsTea, Cd4Te4, CdsTes and CdeTeg, there are two pos- 
sible lowest energy configurations within the energy dif- 
ference of - 0.1 leV, ~ 0.22eV ~ 0.04eV, - 0.05eV and 



~ 0.09eV per atom, respectively. Except for CdgTeg, 
the lowest energy geometry is a planar ring-like struc- 
ture. Cd2Te2 has LE as the rhombic structure and a 
square planar structure as FLM. The LE and FLM be- 
long to D2h and C2v point group repectively. CdsTea 
has triangular structure (-D3/1 point group) as LE config- 
uration and another planar structure {Cs point group) 
as the FLM. No three dimensional (3D) stable geometry 
is obtained for n = 2,3. For Cd4Te4, the square planar 
{D4h point group) geometry is obtained for LE and 3D 
tetrahedral structure {Td point group) as FLM. CdsTes 
is a planar pentagon (Cs point group) in LE and a non- 
symmetric 3D structure as FLM. A transition from pla- 
nar LE to 3D LE occurs for n — 6. This is attributed, as 
explained later, to the tendency of achieving higher co- 
ordination without putting too much strain on the bond 
angles. All the atoms in planar clusters, except for the 
FLM of CdsTes, have coordination 2. 

The geometries obtained for the CdTe clusters are 
quite similar to those of corresponding CdS and CdSe 
clusters upto n = 42i. However, the geometries of 
CdsTes and CdgTeg are very different from the CdS and 
CdSe counterpartsSl. This difference may be attributed 
to the fact that we have included 4d of Cd in the va- 
lence configuration in the present study. The necessity 
of including 4d electrons of the cations in the valence con- 
figuration has been emphasized by Wei and Zunger^^. 

As we move from planar to 3D structure the coordi- 
nation of atoms increases from 2 to 3 and in some cases 
4. The 3D structures of Cd4Te4, CdgTeg and CdeTeg 
can be considered to be made up of smaller cluster units 
namely Cd2Te2 and CdaTca. 

A comparision of Cd„Te„ clusters with Cdrf^i shows 
that three dimensional structures are more favorable in 
pure Cd clusters as compared to CdTe clusters where 
planar structure is favored even for 10 atom cluster. 

Figure 13 shows the variation of average Cd-Te bond 
length with the number of atoms in the clusters. The 
near neighbour (nn) distance for CdTe bulk is 2.81A. As 
evident from Fig. |31 these clusters have an average bond 
length less than the nn distance in bulk. This is due to 
incomplete coordination in these clusters on account of 
very small number of atoms. The bond length of CdTe 
dimer is the smallest. In fact, the smallest Cd-Te bond 
length remains almost constant in these clusters for n = 
3 — 6 (not shown in figure). The average Cd-Te bond 
length shrinks with increasing n in planar structures {n = 
2 — 5) and increases in going from planar to 3D structures. 

Surface studies^^ on hcteropolar covalent and ionic 
semiconductors have identified quantum mechanical 
electron-electron Coulomb repulsion, hybridization ef- 
fects and classical Coulomb attraction between anion and 
cation as competing factors in deciding the stability of 
the surface. The last factor seems to be more important 
with only 3 neighbours for atoms on the surface. Opti- 
mum energy gain is achieved when the distance between 
anion and cation is as small as possible. Thus shorter 
bonds are more ionic in nature. We therefore expect that 



CdTe dimer will have more ionic character as compared 
to other clusters. In Cd2Te2, there exist a covalent bond 
between Cd-Cd which accounts for an increase in average 
Cd-Te bond length. As explained later from the charge 
density plots, there are no covalent bonds between either 
Cd-Cd or Te-Te in other planar structures (n = 3 — 5). 
Therefore, the average Cd-Te bond length is smaller in 
these clusters. In CdgTeg, bonds between similar atoms 
appear and as a result the average Cd-Te bond length 
increases. 

Cd-Cd bond length in Cd2Te2 is smaller (2.871) than 
in Cd dimer (3.47 A). The average Cd-Cd bond length is 
same in CdsTea and Cda. However, it is larger in Cd„Te„ 
as compared to Cd„ for 3<n<6. In contrast, average Te- 
Te bond length is significantly larger in CdTe clusters in 
comparison to pure Te clusters as expected due to lone 
pair repulsion. Also it does not appreciably change with 
n. 

Average Te-Cd-Te bond angle increases from 107.2° 
in Cd2Te2 to nearly 180° in Cd4Te4. The angle further 
opens up in CdsTes. In CdgTee, all the Te-Cd-Te an- 
gles are smaller than linear angle as a consequence of 3D 
structure. Thus we see that smaller clusters have a ten- 
dency to form linear Te-Cd-Te bonds resulting in planar 
structures with coordiantion of 2 and as higher coordina- 
tion number takes over, a transition from planar to three 
dimensional structure results for stable geometry. 

All the clusters are found to be stable against fragmen- 
tation into smaller clusters. 



B. Electronic properties of Cd„Te„ 

Figures 01 and show the variation in the highest oc- 
cupied molecular orbital (HOMO) - lowest unoccupied 
molecular orbital (LUMO) gap (Eg) and binding energy 
(defined in Eq. (^ per Cd-Te unit with the number of 
Cd-Te units in the cluster. 

Eb = [n{E{Cd) + E{Te)} - S(Cd„Te„)]/n (1) 

The corresponding two curves for CdS and CdSe show 
similar trends^^ and the authors claim to infer the stabil- 
ity of clusters from the maximum in these curves. How- 
ever, we do not see any such correlation for the clusters 
studied in the present work. It is seen that CdsTes is 
more stable compared to its neighbours. It may be men- 
tioned that the number of electrons in this system are 90 
which matches with the shell closing number for jellium- 
like systems. 

We have calculated the vertical detachment energy 
(VDE), shown in Fig. |H| and the electron affinity (EA), 
shown in Fig. \7\ using Eqs. |(2Jl and (O by performing 
two self-consistent calculations for different elecronic sys- 
tems without allowing any change in the lowest energy 
geometry of the neutral cluster. 

VDE{CdnTen) = S(Cd„Te„+) - ^(Cd„re„) (2) 



EA{Cdr,Te„) = EiCdnTen') - E{CdnTen) (3) 

In bulk CdTe, the bonding is largely covalent as the 
atomic number and hence the atomic radii for Cd and Te 
are very close. As a result, the Coulombic forces between 
the two species are balanced very well. However due to 
difference in electronegativity of Cd and Te, the bonding 
has some ionic character even in the bulk. 

We show total charge density iso-surface plots for CdTe 
dimer and Cd2Te2 cluster in Figs. El and El respectively, 
revealing covalent bonding. Similar charge density dis- 
tributions are seen for other clusters too. No Cd-Cd or 
Te-Te near neighbour bonds are seen to be formed for 
clusters with n > 2. The charge density plots for CdS 
and CdSe^- show that the bonding has some ionic char- 
acter and is not purely covalent. lonicity in bonding is 
more in CdS than in CdSe because of smaller atomic ra- 
dius of sulphur. lonicity is expected to be still smaller in 
CdTe. 

An analysis of charge inside atomic spheres with cova- 
lent radii shows that the distribution of charge in the Cd 
and Te spheres is almost same for planar structures for 
n — 3 — 5. There is more charge localization in planar 
structures in comparision to CdTe dimer, Cd2Te2 and 
the 3D structures. 

A more detailed information regarding the bonding in 
these clusters can be obtained from the partial charge 
density plots of the molecular orbitals (MOs) in these 
clusters. 

From Figs. [TUl and El it is evident that the HOMO 
consists of purely Te p and in the LUMO Cd s and p 
hybridize to give a sp-hybrid orbital which in turn forms 
a a bond with the Te p orbital. Similar trends are seen 
for the plots for HOMO and LUMO of CdgTca as shown 
in Figs. El and El Figure El also shows the delocal- 
ization of the charge density in the central region of the 
clusters indicating a semi-metallic character in the bond- 
ing. Thus, these systems will favour electron addition as 
is also evident from Fig. [7\ For other clusters the par- 
tial charge density for HOMO and LUMO have similar 
nature. 

As seen from the MO charge distribution of HOMO 
as one goes from n = 2 to n = 6, an enhancement in 
delocalisation results accounting for a steady decrease in 
the VDE as shown in Fig. O 

In Fig. El we show the one-electron energies for all 
the CdTe clusters and the atomic levels of Cd and Te. 
The Cd s and the Te p levels lie in the same energy 
window in the atomic case. Consequently, an intermix- 
ing of these orbitals occurs in the clusters. A fat band 
analysis of the CdTe bulk band structure (not shown 
here) using full potential linear augmented plane wave 
method^iSi reveals similar intermixing of orbitals. Thus, 
the bulk type hybridization of the orbitals is retained in 
clusters. The movement of the the HOMO and LUMO 
in these small clusters is quite interesting. As the clus- 
ter size increases, the HOMO moves downward whereas 



the LUMO moves upward resulting in widening of the 
HOMO-LUMO gap. However, this shift is not uniform. 
The HOMO moves downward slowly as compared to the 
LUMO which moves upward rapidly. But for CdgTeg, 
the HOMO shifts slightly upward and the LUMO shifts 
downward resulting in a smaller HOMO-LUMO gap as a 
consequence of planar to 3D transition in geometry. 



C. Properties of Cd,„Te„ clusters 

CdTe, with many other II- VI semiconductors, in bulk 
has ZB and wurtzite (W) phases possible. There are also 
reports of size induced phase transitions in ZnS^^ and 
CdS^^ nanoparticles. The first coordination for central 
atom is identical in ZB and W phases. If we consider the 
trigonal axis [111] in ZB and [001] in W then three near 
neighbours are identical. 

We have, therefore, performed calculations for four 
larger non-stoichiometric CdTe clusters viz. CdisTcig 
(central atom -I- 3 neighbouring shells), CdigTcig (cen- 
tral atom -I- 4 shells), CdigTcia and CdigTcig. The later 
two geometries are obtained from the former two by inter- 
changing the positions of Cd and Te atoms. These clus- 
ters have been considered as fragments of the bulk and 
hence their initial symmetry is Td as shown in Figs. llSf a) 
and llGr al. 

CdisTcie has a central Cd atom and the side view of 
its relaxed structure is shown in Fig. llSf b'). Figure ITST c') 
shows the side view of relaxed structure of CdieTeia 
which has a central Te atom. During relaxation the cen- 
tral Te atom pushes the neighboring Cd atoms outward. 
Figure llfcil shows the top views of initial geometry for 
Cdi6Tei9(a) and the relaxed geometries for Cdi6Tei9(b) 
and CdigTei6(c). It is found that upon relaxation only 
the Te-rich clusters retain their Td symmetry whereas 
the Cd-rich clusters attain lower symmetry structures. 

We summarize in table |2 some of the structural and 
electronic properites of CdmTe„ clusters. It may be men- 
tioned that the binding energy of CdTe dimer is 0.76 
eV/atom whereas the cohesive energy of CdTe bulk is 
4.9 eV. 



TABLE I: Summary of Structural and Electronic proper- 
ties of relaxed CdmTe„(m7^rz) clusters. PG is the observed 
point group symmetry, BE is binding energy expressed in 
(eV/atom), / is average Cd-Te bond length in A, 9 is the 
Y-X-Y bond angle at the central atom X of the cluster and 
Eg is the HOMO-LUMO gap in eV. 



Clusters 


PG 


BE 


I 


9 


Es 


CdiaTeie 


Td 


2.13 


2.81 


109.47° 


0.195 


CdieTeia 


Csv 


1.95 


2.84 


109.48° 


0.659 


CdieTeig 


Td 


2.14 


2.82 


109.47° 


0.088 


CdigTeie 


Ci 


1.95 


2.81 


91.31° 


0.964 



A comparison of the binding energies of the above clus- 



ters indicates that CdiaTeig and CdigTeig are more sta- 
ble than CdieTeia and CdigTeie- In the former group 
the Te atoms are on the surface and hence these struc- 
tures are energetically favoured as discussed earlier for 
smaller clusters. For the Te-rich clusters, the relaxation 
is not strong enough to break their initial T^ symmetry. 
On the other hand, the Cd-rich clusters, being weakly 
bound, go to lower symmetry structures upon relaxation. 
A comparision of binding energy per atom for CdmTe„ 
with corresponding HgmTe„ clusters^ shows an increase. 
This is attributed to symmetry allowed pd interactions. 
The d levels of Hg are closer than those of Cd to the p 
levels of Te. Therefore pd-repulsion induced bond weak- 
ening is more in HgTe than in CdTe clusters resulting in 
a reduction of binding energy. 

A very interesting feature observed in the Cd-rich clus- 
ters is Jahn-Teller distortion. The initial HOMO-LUMO 
gap for Cdi6Tei3 and CdigTeie is nearly zero. But upon 
relaxation, they attain an appreciable HOMO-LUMO 
gap. This is a clear indication that the energy levels 
which are nearly degenerate in the initial structure have 
moved apart during relaxation. Thus, relaxation has de- 
stroyed symmetry resulting in the lifting of degeneracy. 
On the other hand for the Te-rich clusters, the HOMO- 
LUMO are not degenerate (though close to each other) 
and thus there is no Jahn-Teller distortion. Similar trend 
has also been observed in non-stoichiomctric HgTe clus- 
ters*. The Hg-rich clusters upon relaxation lose their 
initial T^ symmetry. As a result, they attain a definite 
bandgap and a semi-metal to semiconductor transition 
is observed. Dalpian et al. have shown that symmetry 
is an important parameter in deciding the properties of 
nanoclusters^*. 



resemble that of the bulk phase. The Cd atoms tend to 
form a core group surrounded by the chalcogenide atoms 
because of the Coulomb repulsion of the lone pair of p- 
electrons on Te. For Cd„Te„, (l<n<5), the planar geom- 
etry is the lowest energy configuration as Te-Cd-Te prefer 
to form linear bond. But as the number of atoms in- 
creases further, higher coordination takes over and three 
dimensional structures result. Bonding in CdTe clusters 
is mostly covalent in nature with partial ionic character. 
The degree of ionicity changes with n. Hybridization of 
orbitals in these clusters is observed to be same as that 
in the bulk. 

The relaxed structures of non-stoichiometric clusters, 
Cd„iTe„ with (m,,n = 13, 16, 19 and m ^ n), were ob- 
tained using ab-initio DFT with the wavefunction pro- 
jected to real space. The Te-rich clusters with Te atoms 
on the surface, are found to be more stable than the Cd- 
rich clusters. The symmetry of the cluster plays a vital 
role in determining its HOMO-LUMO gap. It was found 
that the Cd-rich cluster lose their initial T^d symmetry 
and attain a large HOMO-LUMO gap as compared to 
the Te-rich clusters. 



ACKNOWLEDGMENT 



IV. CONCLUSION 

The ground state geometries of Cd„Te„ (l<n<6) were 
calculated using ab-initio DFT. The structures do not 



We acknowledge financial support from the Depart- 
ment of Science & Technology, Governement of India. 
We also thank lUCAA, Pune and Centre for Modeling 
and Simulation, University of Pune for use of their com- 
putational facilities. 



* Electronic address: |skb@physics . unipune . ernet . in| 

t Electronic address: "anjali@phy sics. unipune. ernet. in| 

^ Corresponding author. 
^ V. L. Colvin, M. C. Schlamp and A. P. Alivisatos, Nature 

(London) 370, 354 (1994); A. P. Alivisatos, J. Phys. Chem. 

100, 13226 (1996). 
^ S. Ogot, J. R. Chelikowsky, S. G. Louie, Phys. Rev. Lett. 

79 1770 (1997); L Vasiliev, S. Ogot and J. R. Chelikowsky, 

Phys. Rev. B 65, 115416 (2002) and references therein. 
^ A. Shvartsberg (private communication). 

* M. Gao, C. Lesser, S. Kirstein, H. M/"ohwaId, A.L. Ro- 
gach and H. Weller, J. Appl. Phys. 87, 2297 (2000). 

^ Li J ei al, Chem. Commn. (Camb.) 8, 982 (2004). 

® N. W. Ashcroft and N. David Mermin, Solid State Physics, 
(Holt-Saunders, Philadelphia, 1981). 

"^ M. L. Cohen and J. R. Chelikowsky in Electronic Struc- 
ture and Optical Properties of Semiconductors edited by 



M. Cardona, (Springer- Verlag, Berlin, 1988). 
* X. Q. Wang, S. J. Clark and R. A. Abram, Phys. Rev. B 

70, 235328 (2004). 
® V. S. Gurin, Int. J. Q. Chem. 71, 337 (1999); J. -O. Joswig 
et al, J. Phys. Chem. B 104, 2617 (2000); A. J. WiUiamson 
and A. Zunger, Phys. Rev. B 61, 1978 (2000); J. M. Matx- 
ain, J. E. Fowler and J. M. Ugalde, Phys. Rev. A 61, 
053201 (2000); J. M. Matxain, J. M. Mercero, J. E. Fowler 
and J. M. Ugalde, Phys. Rev. A 64, 053201 (2001). 

1° J. Muilu and T. A. Pakkanen, Sur. Sci. 364, 439 (1996). 

^^ E. Spano, S. Hamad and C. R. Catlow, J. Phys. Chem. 13, 
10337 (2003). 

^^ V. Albe, C. Jouanin and D. Bertho, J. Crys. Gr. 184/185, 
388 (1998); S. Pokrant and K. B. Whaley Eur. Phys. J. D 
6, 255 (1999). 

^^ P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964); 
W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 



" p. E. Blochl, Phys. Rev. B 50, 17953 (1994); G. Krese and 
J. Joubert, Phys. Rev. B 59, 1759 (1999). 

^^ G.Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 
(1996); Comput. Mater. Sci. 6, 15 (1996). 

^^ A. Pasquarello, K. Laasonen, R. Car, C. Lee and D. Van- 
derbilt, Phys. Rev. Lett. 69 1982 (1992); K. Laasonen, A. 
Pasquarello, R. Car, C. Lee and D. Vanderbilt, Phys. Rev. 
B 47, 10142 (1993). 

^^ O. Gunnarsson and R. Jones in Local Density Approxi- 
mations in Quantum Chemistry and Solid State Theory, 
edited by J. P. Dahl and J. Avery, (Plenum, NewYourk, 
1983), pp. 229. 

^® J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. 
77, 3865 (1996). 

^^ J. P. Perdew Electronic Structure of Solids '91, edited by 
P. Ziesche and H.Eschrig (Akademie Verlag, Berlin, 1991), 
pp. 11; J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 
(1992). 



2° N. Kosugi, J. Comput. Phys. 55, 426 (1984). 

^^ M. C. Troparevsky and J. R. Chelikowsky, J. Chem. Phys. 
114, 943 (2001). 

^^ S. -H. Wei and A. Zungcr, Phys. Rev. B 37, 8958 (1988); 
J. Cry. Gr. 86, 1 (1988). 

2^ J. Zhao, Phys. Rev. A 64, 043204 (2001). 

^^ D. Vogel, P. Kriigcr and J. Pollmann, Sur. Sci. 402-404, 
774 (1998). 

^^ P. Blaha, K. Schwarz, P. Sorantin and S.B. Trickcy, Comp. 
Phys. Commun. 59, 399 (1990). 

2^ J. Nanda and D. D. Sarma, J. Appl. Phys. 90, 2504 (2001). 

^■^ T. Nanba, Y. Nodake, M. Muneyasu, G. P. Wilhams and S. 
Hayashi, J. Phys. See. Jpn. 66, 1526 (1997); R. Banerjee, 
R. Jayakrishnan and P. Ayyub, J. Phys. : Cond. Matt. 12, 
10647 (2000). 

^* G. M. Dalpian, M. L. Tiago, M. L. Puerto and J. R. Che- 
likowsky, Nano Lett. 6, 501 (2006). 



■^ 



^ 






FIG. 2: Geometries of the first local minima of CdnTen for 
(1 <n< 6). 






c 

c 
o 
x^ 
m 
H 
-a 
U 






2.85 



2.75 



2.65 



2.55 




12 3 4 5 6 

no. of Cd-Te units 

FIG. 3: Average Cd-Te bond length in A vs number of Cd-Te 
units. 
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FIG. 4: HOMO-LUMO gap in eV vs number of Cd-Te units. 
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FIG. 5: Binding energy in eV vs number of Cd-To units. 
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FIG. 6: Vertical detachment energy (eV) vs number of Cd-Te 
units. 
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FIG. 7: Electron affinity in eV vs number of Cd-Te units. 







FIG. 8: Charge density iso-surface for CdTe dimer. 




FIG. 9: Charge density iso-surface for Cd2Te2. The grey 
atoms are Cd and black atoms are Te. 
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FIG. 10: Iso-surface of partial charge density of HOMO for 
CdTe dimer. 
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FIG. 11: Iso-surface of partial charge density of LUMO for 
CdTe dimer. 
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FIG. 12: Iso-surface of partial charge density of HOMO for 
CdsTea. 




FIG. 13: Iso-surface of partial charge density of LUMO for 
CdsTea. 
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FIG. 14: Energy level diagram of CdTe clusters. 



jjj''^§^*3 




Y v*- •■ ■*, w 




(a) CdisTeie : initial 



(b) CdisTeie '■ relaxed (c) CdieTeis : relaxed 

FIG. 15: Side view of the initial and relaxed geometries. 






(a) CdiaTeig : initial 



(b) CdieTeiQ : relaxed (c) CdigTeie : relaxed 

FIG. 16: Top view of the initial and relaxed geometries. 



